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We demonstrate the emission of highly indistinguishable photons from a quasi-resonantly pumped 
coupled quantum dot-microcavity system operating in the regime of cavity quantum electrodynam¬ 
ics. Changing the sample temperature allows us to vary the quantum dot-cavity detuning, and on 
spectral resonance we observe a three-fold improvement in the Hong-Ou-Mandel interference visibil¬ 
ity, reaching values in excess of 80%. Our measurements off-resonance allow us to investigate varying 
Purcell enhancements, and to probe the dephasing environment at different temperatures and en¬ 
ergy scales. By comparison with our microscopic model, we are able to identify pure-dephasing and 
not time-jitter as the dominating source of imperfections in our system. 


I. INTRODUCTION 

Single indistinguishable photons are key to appli¬ 
cations in quantum networks [T] , linear optical quan¬ 
tum computing mm and quantum teleportation HH^- 
One of the most promising platforms for single pho¬ 
ton sources are solid-state quantum dots (QDs) |SlfIUj. 
Compared to alternative platforms, such as cold atoms 
or trapped ions, single QDs offer several advantages: 
they can be driven electrically, which is of crucial im¬ 
portance for compact future applications hhe], and 
in principle can be integrated in complex photonic en¬ 
vironments and architectures, such as on-chip quan¬ 
tum optical networks mm]. When embedded in a 
bulk semiconductor, however, QDs suffer from poor 
photon extraction efficiencies, since only a minor frac¬ 
tion of the photons can leave the high refractive index 
material. This problem can be mitigated by integrat¬ 
ing QDs into optical microcavities rnrnm or pho¬ 
tonic waveguides fTTH-P 11 . which can enhance extrac¬ 
tion efficiencies to values beyond 50%. 

In addition to increased extraction efficiencies, ex¬ 
ploiting cavity quantum electrodynamics (cQED) ef¬ 
fects in QD-based sources can have a positive effect 
on the interference properties (and hence the indis- 
tinguishability) of the emitted photon wave packets. 
Ideally, the wave packets emitted by an indistinguish¬ 
able photon source are Fourier-limited, with a recom¬ 
bination time Ti , and temporal extension of the wave 
packet given by T 2 = 2Ti ( 22 ] ■ If additional dephasing 
channels with a characteristic time l/y exist, such as 
coupling to phonons or spectral diffusion, the coher¬ 
ence time is reduced according to ^ = ^jr- + 7, which 
consequently leads to a reduction of the two photon 
interference visibility. It was theoretically shown that 
pure dephasing strongly affects the detuning depen¬ 
dence of the relative strength of the cavity and QD- 
emission peaks [221 123 ] . In the regime of cQED, the 
lifetime of the QD excitons can be manipulated via 
the photonic density of states in the cavity (the Pur¬ 


cell effect). If the timing of emission events is precisely 
known, and 7 is constant, shortening of the emitter 
lifetime T) via the Purcell effect leads to an improved 
interference visibility as the condition T 2 = 2T\ can 
be approximately restored This simple 

picture, however, is known to breakdown if there are 
uncertainties in the timing of emission events (time- 
jitters) [26H28I . or if the dephasing environment gives 
rise to more than a simple constant pure-dephasing 
rate, as is known to be the case for phonons [251432] . 
As such, with the aim of designing improved single in¬ 
distinguishable photon sources, it is crucially impor¬ 
tant to first establish the magnitude of time-jitters 
and the nature of any dephasing environments. 

In this work, we exploit a microcavity with a high 
Purcell factor and weak non-resonant contributions of 
spectator QDs to probe the interference properties of 
photons emitted from a single QD as a function of the 
QD-cavity detuning. In contrast to previous studies, 
where non-resonant coupling to spectator QDs [33] or 
strong temperature induced dephasing [25] dominated 
the experiments, we observe a strong improvement of 
the two-photon visibility on resonance, which exceeds 
a factor of 3 compared to the off-resonant case. We 
extend the theoretical model of Ref. [25] to derive an 
expression for the Hong-Ou-Mandel dip including the 
effects of both tinre-jitter and pure-dephasing on- and 
off-resonance. This allows us to reject timing-jitter, 
and definitively attribute sources of pure-dephasing 
as the dominant factor limiting the indistinguishabil- 
ity of our photons. Furthermore, we show that the 
degree of symmetry we observe for positive and nega¬ 
tive detuning suggests pure-dephasing caused by both 
phonon coupling and spectral diffusion. 


II. QUANTUM DOT-CAVITY SYSTEM 

The device under investigation comprises a QD em¬ 
bedded in a micropillar cavity with a quality factor 
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FIG. 1: a) Temperature dependent intensity map of the 
QD-cavity system. The QD reaches spectral resonance 
with the cavity mode at T = 17.5 K. b) QD lifetime as a 
function of QD-cavity detuning. The fit is a Lorentzian 
profile where the linewidth is fixed to the cavity linewidth. 
A Purcell enhancement of Fp = 7.8 ± 2.3 is extracted. 


of Q = 3200. The layer structure consists of 25 (30) 
alternating j^-thick GaAs/AlAs mirror pairs which 
form the upper (lower) distributed Bragg reflector 
(DBR). The cavity region is composed of six alter¬ 
nating GaAs/AlAs layers with decreasing (lower part) 
and increasing (upper part) thickness. A single layer 
of partially capped and annealed InAs QDs is inte¬ 
grated in the central layer of the tapered segment, i.e. 
in the vertical maximum of the optical field [M|. Mi¬ 
cropillars with varying diameters were etched into the 
wafer (the pillar under investigation has a diameter of 
1050 nm) to provide zero dimensional mode confine¬ 
ment. As a result of the Bloch mode engineering [53] , 
our micropillars support optical resonances with com¬ 
parably large quality factors down to the sub-micron 
diameter range, which yields the possibility to signif¬ 
icantly increase the Purcell factor in such microcavi¬ 
ties compared to conventional DBR resonators based 
on A-thick cavity spacers. The sample was placed in¬ 
side an optical cryostat, and the QD was excited via 
a picosecond-pulsed Ti:sapphire laser with a repeti¬ 
tion frequency of 82 MHz (pulse separation 12.2 ns). 
The laser beam was coupled into the optical path via 
a polarizing beam splitter, which also suppresses the 
scattered laser light from the detection path of the 
setup. Further filtering was implemented by a long- 
pass filter in front of the monochromator. After spec¬ 
tral filtering, the emitted photons were coupled into a 
polarization maintaining single mode fibre followed by 
a fibre coupled Mach-Zehnder-Interferometer (MZI) 
with a variable fibre-coupled time delay in one arm to 
measure the two-photon-interference in a Hong-Ou 
Mandel (HOM) setup. The second beamsplitter of the 
MZI can be removed to directly measure the autocor¬ 
relation function of the signal. 

Fig. 0 (a) shows the temperature dependent micro¬ 
photoluminescence (/r-PL) map of the investigated 
QD-cavity system, which was recorded under non¬ 
resonant excitation conditions. The QD emission line, 
which we attribute to the neutral exciton, can be 
tuned through the cavity mode by changing the sam¬ 
ple temperature. Spectral resonance with the funda¬ 
mental cavity mode is achieved at T = 17.5 K. Due 


to the Purcell enhancement, the integrated intensity 
of the QD increases by a factor of more than three 
when the QD and cavity are tuned into resonance. 
In order to directly and accurately extract the Pur¬ 
cell factor of our coupled system, we measured the 
exciton lifetime via time-resolved /i-PL as a function 
of the QD-cavity detuning (see Appendix A2). As 
seen in Fig. [l] (b), we observe a strong decrease of 
the lifetime when the QD is tuned into resonance 
as a result of the Purcell effect. The Purcell factor 
F P = T Ti{a^o) ~ 1 05. * s extracted by fitting the 
data with a Lorentzian profile (the width being fixed 
to the cavity linewidth k = 0.42 meV), and we find a 
value as high as Fp = 7.8 ±2.3 as a result of the small 
mode volume of our microcavity. 


We now study the single photon emission proper¬ 
ties of our system, which is particularly important 
on-resonance, where the single photon characteristics 
can be deteriorated by non-resonant contributions to 
the cavity from spectator QDs, or luminescence from 
the background continuum funnelled into the cavity 
mode [36]. The second-order photon-autocorrelation 
was probed under quasi-resonant excitation condi¬ 
tions, with a laser tuned 32 meV to the high energy 
side of the single exciton emission feature, with a (be¬ 
low saturation) power of 311 fj, W. The on-resonance 
(T = 17.5 K) autocorrelation histogram is shown in 
Fig. [2] (a). The strongly suppressed peak around r = 0 
is a clear signature of single photon emission. We 
extract the g^ 2 \r = 0) value by dividing the area 
of the central peak by the average area of all the 
side peaks, leading to g^ 2 \ 0) = 0.058 ± 0.006, re¬ 
flecting the high purity of our cavity enhanced sin¬ 
gle photon source. Off-resonance we find a minimum 
value of ^ 2 >(0) = (0.036 ± 0.005) at A = -0.28 meV 
(T = 6.4 K). For increasing temperatures, we note 
a modest increase up to g (2 ^(0) = (0.11 ± 0.01) for 
A = 0.61 meV (T = 25.5 K). This value is still close 
to perfect single photon emission, and we attribute 
the slight rise to a lowered signal to background ra¬ 
tio between QD and cavity emission. We note that 
no deterioration of the g^ 2 \0 ) value can be observed 
on spectral resonance, which suggests only very weak 
contributions from spectator QDs to the cavity signal. 




-0.4 -0.2 0.0 0.2 0.4 0.6 0.8 

Detuning (meV) 


FIG. 2: a) Autocorrelation histogram on QD-cavity res¬ 
onance, from which we extract </ 2) (0) = (0.058 ± 0.006). 
b) Autocorrelation function as a function of QD-cavity 
detuning (the inset shows a zoom-in of the y-axis). 
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FIG. 3: a) Histogram of the two-photon-interference for 
zero time delay between the two arms of the interferom¬ 
eter, from which we extract 5 hom( 0) = 0.17 ± 0.02. b) 
Measured two-photon-interference versus the time delay. 
The measurement shows the clear HOM-dip (black dots). 
Due to a non-vanishing two-photon probability, the data 
points go slightly above 0.5 for large td ■ These are cor¬ 
rected by subtracting half the measured </ 2 l(0) value on 
resonance shown in Fig. [2] (red markers). These corrected 
data were fitted to extract a visibility of v = (83 ± 5)%. 


III. PHOTON INDISTINGUISHABILITY 


We now assess the indistinguishable nature of the 
emitted photons, which we probe in the HOM inter¬ 
ferometer under the same pulsed quasi-resonant ex¬ 
citation conditions. The second order correlation his¬ 
togram for zero time delay in the MZI for the resonant 
case is shown in Fig. [3] (a). Strong suppression of the 
central correlation peak directly reflects a strong de¬ 
gree of photon indistinguishability. The black markers 
in Fig. [3] (b) are obtained by dividing the area of the 
peak centred around r = 0 by that centred around 
t = 24 ns for various time delays td, and we observe 
a clear HOM-dip. For large time delays t b the corre¬ 
lation values slightly exceed 0.5 as a result of the finite 
two photon emission probability, as seen in Fig. [2j We 
correct the interference data by subtracting half the 
corresponding experimentally extracted on-resonance 
value of (0 ) = (0.058 ± 0.006) [red markers] (see 
Appendix B 2). We then fit our data to the func¬ 
tion ShomHd) = 0.5(1 - i/exp[—|TD|/Ti]), where we 
set Ti = 67 ps (see Fig. [l] (b)), and we find a visibility 
of v = 83 ± 5%. This high value is a direct conse¬ 
quence of the large Purcell factor in our high quality 
QD-cavity system. 


To further analyse our experimental data, and in 
particular, to determine the relative influences of 
time-jitter and pure-dephasing on the indistinguisha¬ 
bility of the emitted photons, we extend the theory of 
Ref. [28] to derive an expression for the TPI as a func¬ 
tion of both time delay tjj and detuning. Dephasing 
caused by coupling to phonons is known to affect the 
two-photon interference (TPI) properties of the emis¬ 
sion from a QD-cavity system in a highly non-trivial 
way, giving rise, for example, to pronounced asym¬ 
metries for positive and negative QD-cavity detun¬ 
ings [28H30] . We find, however, that nearly all features 
seen in our data can be well reproduced by a model 
assuming a simple constant pure-dephasing rate. We 


present this simplified model first, and then go on to 
show that by including phonons in a rigorous manner 
at a Hamiltonian level, the behaviour off-resonance 
allows us to approximately determine the relative in¬ 
fluence of phonons as compared to other sources of 
dephasing. 

We model the QD as a three-level-system, and con¬ 
sider the vacuum and single photon Fock states of the 
cavity. Provided the QD-cavity coupling strength is 
sufficiently weak, and/or the cavity decay rate is suf¬ 
ficiently large, the cavity degrees of freedom can be 
adiabatically eliminated from the equations of motion 
for the QD-cavity system [28] . The result is a master 
equation of the form (see Appendix |B 3[ ) 

f t =-i[A\E)(E\,p]+ 

{L r (\G)(E\) + L 2l (\E)(E\) + L a (\E)(P\))p, (1) 

where the states | E) = |e,n = 0), |G) = |g,n=l) 
\P) = |p, n = 0) represent the QD in ground (g), single 
exciton state (e), or pump-level (p), with the cavity 
containing zero or one excitations. The QD-cavity 
detuning is A, while 7 is the pure-dephasing rate, and 
a is the rate at which the pump-level decays into the 
single exciton state, with T a = 1/a determining the 
magnitude of the time-jitter (i.e. T a = 0 represents 
the ideal case in which there is no time-jitter). The 
Purcell enhanced spontaneous emission rate is 

r = Tr 1 = r B + 2 5 2 -^^, ( 2 ) 

with r B the background decay rate, g the QD-cavity 
coupling strength, and 7 t 0 t = 7 + ^(k + with k 
the cavity decay rate. The validity of Eq. ([!]) relies on 
the condition 7 to t A,T,g, which is satisfied in all 
our experiments. 

Eq. 0 can then be used to derive an expression 
for the normalised coincidence events in the TPI mea¬ 
surements (for details see Appendix |B 3[ ). The second 
order correlation function for the HOM interference 
measurements is found to read 

ShomM = 

where the detuning dependence enters through V [see 
Eq. 0], and u = (r/(r+ 27 ))(a/(r+a')) is the visibil¬ 
ity. We note that while the expression for v has been 
derived before [23, to our knowledge Eq. 0 repre¬ 
sents the first time the full behaviour of the HOM- 
dip for nonzero values of including time-jitter and 
pure-dephasing has been presented. This model pro¬ 
vides us with simple analytical expressions with which 
we can fit the experimental TPI data. Crucially, it al¬ 
lows us to explore how a given set of parameters simul¬ 
taneously affects the HOM-dip and the TPI visibility 
as the QD and cavity are moved off resonance. 

In Fig. [i] (a) we show again the HOM-dip, while in 
Fig. [i] (b) we show the depth of HOM-dip as a func¬ 
tion of detuning. We see a pronounced rise of the 
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FIG. 4: a) HOM-dip and HOM-dip depth as a function 
of detuning (b). The dashed blue curves correspond to a 
parameter set for which time-jitter dominates, which we 
find to be inconsistent with the data off resonance shown in 
(b). The dotted green curves correspond to a parameters 
set for which pure-dephasing dominates, which is able to 
consistently reproduce all data. The solid orange curves 
again correspond to a parameter set dominated by pure- 
dephasing, but where 40% of the dephasing on-resonance 
is caused by coupling to phonons. 


HOM-dip as the QD is brought out of resonance, cor¬ 
responding to visibilities on- and off-resonance which 
differ by more than a factor of 3. The dashed blue 
curves in Fi g. [4] show a fit to Eq. ([3]) where only the 
data in Fig. [4[ (a) (the HOM-dip) is considered. Fit¬ 
ting parameters T a = 1/a = 21.5 ps and 7 = 0 /xeV 
are found, corresponding to a regime where time-jitter 
dominates. These parameters are able to reproduce 
the HOM-dip well, but fail to describe the data off 
resonance. Indeed, we find that if we simultaneously 
fit all the data shown in Fig. [4j we find parameters 
T a = 1/a = 3.2 ps and 7 = 1.6 geV, corresponding to 
a regime for which pure-dephasing dominates. These 
parameters are shown by the dotted green curves, and 
much better agreement is found. These two fits show 
that while both time-jitter and pure-dephasing affect 
the shape of the HOM dip in a similar way, the re¬ 
duction in the visibility seen off-resonance can only 
be explained by a system in which pure-dephasing 
dominates. As the QD and cavity are moved off res¬ 
onance, the Purcell effect weakens (see Fig. |T]) and 
r = 1/Ti decreases. For a source dominated by pure- 
dephasing, the visibility is given by v ~ T/(r + 2y), 
and a reduction in T causes a reduction in v. For a 
source dominated by time-jitter, the visibility instead 
follows v ss a/(r + a), and a reduction in F increases 
v. We stress that which of two regimes is relevant 
for a particular system has important consequences 
for how experimental improvements will translate to 
improvements in photon indistinguishabilities. In the 
present case, since pure-dephasing dominates, a com¬ 
plete elimination of time-jitters (achieved, for exam¬ 
ple, via strictly resonant excitation conditions), will 
lead to only a modest 4% increase in the visibility, 
while an elimination of sources of pure-dephasing will 
lead to an increase of 20% up to v = 95%. 


IV. DISCUSSION 


The low value of T a = 3.2 ps implies that our quasi¬ 
resonant excitation scheme leads to a very fast relax¬ 
ation to the desired single exciton state. This is also 
supported by the laser detuning we use (32 meV), 
which corresponds to the energy of a longitudinal op¬ 
tical phonon, known to relax on this timescale |37| . 
We attribute pure-dephasing in our sample as caused 
by exciton-phonon coupling and spectral fluctuation 
of the QD energy levels on a timescale shorter than 
the pulse separation of 12.2 ns. The constant pure- 
dephasing rate used in our theory is expected to well 
approximate the spectral fluctuations, but the influ¬ 
ence of phonons is known to give rise to more compli¬ 
cated behavior [281431] . In particular, differing phonon 
absorption and emission rates at low temperatures are 
expected to lead to asymmetries for positive and neg¬ 
ative detuning ;23|. By including phonons using a 
time-convolutionless master equ ation technique (see 
e.g. Ref. [28] or Appendix B5), we find that these 


asymmetries can improve our fits. The solid orange 
curves in Fig. [4] show the predictions of a parame¬ 
ter set similar to that of the dotted green curve, but 
where we have included phonons with a strength cor¬ 
responding to approximately 40% of the total pure- 
dephasing on-resonance [42], and it can be seen that 
the phonon contribution improves the fits to the data. 
We note, however, that when increasing the phonon 
contribution yet further, the fits become worse as the 
asymmetry becomes too strong. The relatively strong 
symmetry seen in Fig. [4] (b) therefore leads us to con¬ 
clude that both phonons, and additional sources of 
constant pure-dephasing (such as a spectral diffusion) 
are present in our system. 

In conclusion, we have demonstrated the feasibil¬ 
ity of our novel cavity design to enhance the emis¬ 
sion of indistinguishable single photons generated in 
epitaxially grown InAs-QDs by a quasi-resonant ex¬ 
citation scheme with a TPI-visibility as high as v = 
(83 ± 5)%, and a two-photon emission probability as 
low as g^ 2 \ 0) = (0.036 ± 0.005). We studied the in¬ 
fluence of the QD-cavity detuning on both the two- 
photon-probability and the degree of indistinguisha- 
bility of the emitted photons. The TPI measure¬ 
ments are explained by our new theory which takes the 
QD-cavity-detuning, time-jitter and pure-dephasing 
into account, and which identifies sources of pure- 
dephasing as the ultimate factor limiting the indis- 
tinguishably of emitted photons. 
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Appendix A: Experimental Methods 


1. Coherence Measurements 


In addition to the Hong-Ou-Mandel (HOM) inter¬ 
ference measurements, the coherence of the emitted 
photons was measured using a free-beam unbalanced 
Michelson interferometer. One mirror is mounted on 
a 300 mm long linear stage which defines the path 
length difference between both optical arms, and using 
an additional implemented piezo crystal at one mirror, 
the contrast of the emitted photons is measured as a 
function of the path length difference. The measure¬ 
ments (black data points) are shown in Fig. A.l (a) 


for a QD in spectral resonance with the cavity mode. 
Fitting these data points to a Gaussian function of 
the form A + B * exp[— (7 t/2)(t/T2 ) 2 ] we extract a 
coherence time of T 2 = (93 ± 3) ps. This value is 
slightly lower than the coherence time extracted from 
the HOM-dip in Fig. [3j for which T 2 = 111 ps. We 
attribute this slight discrepancy to a long term spec¬ 
tral jitter which affects the QD emission energy on 
timescales which are longer than the pulse separation. 
In the HOM-measurement, only subsequently emitted 
photons separated by 12.2 ns (the laser pulse separa¬ 
tion) contribute to the measured indistinguishability, 
and hence the inferred coherence time of T 2 = 111 ps. 
The HOM measurements therefore include an effective 
time filter. In contrast, the measurements made using 
the Michelson interferometer are time-integrated, and 
as such long-term drifts and spectral diffusion result 
in a deterioration of the extracted T 2 value mmm- 


2. Lifetime Measurements 


In order to measure the lifetime of the QD emis¬ 
sion we couple the spectrally filtered photons into 
a single mode fibre attached to an avalanche photo 
diode (APD) with resolution ~ 40 ps. Fig. A.l (b) 
shows two representative time-resolved measurements 
of the QD emission under quasi-resonant excitation. 
The blue round data points correspond to spectral 
resonance between QD and fundamental cavity mode 
(A = 0 meV), while the red square data points 
correspond to a detuning of A = 0.61 meV. The 
measurements (time window 100 ns) each contain 
six complete decay curves similar to those shown in 
Fig. A.l (b), which we fit to a biexponential decay 
function. The shorter time constant represents the 
lifetime of the briglit-exciton, while the longer origi- 


FIG. A.l: a) Contrast of the emission measured using a 
free-beam unbalanced Michelson interferometer as a func¬ 
tion of path length difference of both arms. Fitting the 
data points with a Gaussian distribution (red solid curve) 
we extract a coherence time of T 2 = (93 ± 3) ps. b) Time- 
resolved /i-PL measurements. The blue dotted data points 
correspond to spectral resonance between QD and funda¬ 
mental cavity mode (A = 0 meV), while the red square 
points were taken for a detuning of A = 0.61 meV. The 
solid curves correspond to fits to a biexponential function, 
from which we extract T\ = (67 ± 8) ps on resonance and 
Ti = (306 ± 13) ps off resonance. 


nates from a dark exciton effect. For the decay curves 
in Fig. A.l (b) we find T) = (67 ± 8) ps on resonance 
and Ti = (306 ± 13) ps for A = 0.61 meV. 


Appendix B: Two-Photon Interference Theory 

Here we provide the necessary background for the 
theoretical analysis of the data presented in the main 
text. 


1. Hanbury Brown and Twiss Measurements 

We first consider the Hanbury Brown and Twiss 
(HBT) experimental setup used to measure the two- 
photon emission probability of our source. Emission 
from the source is incident upon a 50/50 beam splitter, 
and two detectors are placed equidistantly on the two 
output arms. We label t\ the time of the detection 
event at detector 1, and t 2 that of detector 2. The 
probability of detecting a photon at detector 1 at t±, 
and at detector 2 at t 2 is proportional to the second 
order field correlation function 

G 1 ' 2 ' 1 (t\,t 2 ) = {b\{ti)b\(t 2 )b 2 (t 2 )b\(t \)), (Bl) 

where b\ (t) is the creation operator for the mode prop¬ 
agating to detector 1 in the Heisenberg picture, and 
similarly for b 2 (t). We relate these modes to those 
on the input arms, described by creation operators a\ 
and a\i using the unitary mode transformation [26 1 



where td is the delay introduced between arrival 
times at the beam-splitter. For the HBT measure¬ 
ment, there is no input in arm 1, and we simply have 
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G (2) (U,f 2 ) -> G < £l T (t 1 ,t 2 ) with 


2. Hong-Ou-Mandel Experiment 


G ( al T (t 1 ,t 2 ) = ^(cl 1 (ti)a 1 (t 2 )a(t 2 )a(ti)) , (B3) 

where the subscripts on the operators have been 
dropped since they are all equal. 

To measure the two-photon emission probability, 
<? (2) (0), we integrate Eq. ( |B3| ) over all t\ and t 2 , and 
divide this area by an adj acen t peak. The adjacent 
peaks correspond to Eq. (B3|, but where t\ and t 2 


differ sufficiently that mode operators at these times 
are completely uncorrelated. This gives the uncor¬ 
related coincidence probability in the HBT measure¬ 
ment G f HB T , UC (* 1 » *a) = (l/4)5 (2) (ii:i2) with 

G {2) {ti,t 2 ) = (a t (ti)a(ti))(a t (t 2 )a(t 2 ))- (B4) 

The normalised autocorrelation function is then de¬ 
fined as 

^( 2 )^ _ f-oo dil f-oo dt 2 GHB T (t 1 ,t 2 ) 


Loo dtl /_OO ^G^ T UN (tl,t 2 ) 

IZ dti fZ dt 2 (a\t\)a)(t 2 )a(t 2 )a{t{)) 


(B5) 


IZ dt i fZ dt 2 (a1(t 1 )a(t 1 ))(a1(t 2 )a(t 2 )) 

which is equal to zero for {a* {t\)(A {t 2 )a{t 2 )a(ti)) = 0. 


We now consider the Hong-Ou Mandel (HOM) ex¬ 
perimental setup used to measure the indistinguish¬ 
able nature of the emitted photons. Two emission 
events are incident on a 50/50 beam-splitter, with 
a delay tjj introduced into input arm one. The un¬ 
normalised probability of a coincidence event is again 
given by Eq. (Bl), and the beam-splitter is described 


by Eq. (B2). Upon combining these equations we find 


16 terms. These can be simplified by assuming that 
modes 1 and 2 are identical but statistically indepen¬ 
dent, which allows us to write ( A 1 A 2 ) = (A 1 )(A 2 ), 
where A\ is any product of mode operators pertain¬ 
ing to mode 1, and similarly for A 2 . We then find eight 
terms linear in (ai) and ( a 2 ). For an electromagnetic 
field state of the form Yl n a ™ l n )( n l’ with | n) a Fock 
state, expectation values linear in the ladder operators 
are zero, and we neglect these terms. This leaves sec¬ 
ond and fourth order terms. The second order terms 
involve expectation values of the form (a) {t\)a} {t 2 )), 
which also give zero for electromagnetic fields as dis¬ 
cussed above. The remaining terms give 


^2) T~d) — GhbT^ 1 t dA2 T £>) + G bbt (U, i 2 ) 

+ | (g {2 \t 1 - T D ,t 2 ) + g { 2 ) (h,t 2 - T D ) - 2Re G (1) (G - r D ,t 2 - t d )G w ( t 2 ,ti) 


-.( 2 ) 


( 2 ) 


(B6) 


where G (1) (fi,f 2 ) = (a t (ti)a(t 2 )) is the unnormalised 
first order correlation function. 

To normalise this quantity we again consider the 
scenario in which t\ and t 2 are sufficiently separated 
that mode operators evaluated at these two times 
are uncorrelated. In doing so we find the uncor¬ 
related coincidence probability for the HOM setup 


G 


( 2 ) 
HOM.UC 


(ti,t 2 ) = \{g ( ' 2) (t 1 ,t2) + - r D ,t 2 ) + 

G^HhA 2 - t d ) + g {2) (ti - T D ,t 2 - t d )). Since we 
integrate over all t\ and t 2 , the appearances of td can 
be neglected, i.e. we have 


/ OO nOO 

dt\ / dt 2 GyQ MUC (ti, t 2 ) = 

-OO J —OO 

/ OO pOO 

dt 1 / dt 2 g (2 \t 1 ,t2) (B7) 

-oo J —oo 


by a simple change of variables. An identical argu- 
ment can be made for t he td appearing in Gjj BT (fi — 
t d , t 2 — td ) in Eq. (B61. In the HOM setup we there¬ 


fore measure the normalised quantity 

IZ d G IZ dt 2G ( Z M (ti, t 2l T D ) 


9hom( t d) — 


fZ^fZ^G^ 


HOM,UN 


(tl,t 2 ) 


= \Z ] {V) + \(} -C(t d )) 


(B8) 


where we have defined the strictly-two-photon coales¬ 
cence probability 

C(t d ) = 

IZ dtl fZ dt 2 Re l G(1) ( tl - t dA2 - UD)G (1) (t 2 ,t i)] 


JZ dt ZZ dt2 S {2) (h,t2) 


(B9) 


which for td = 0 becomes the visibility v = G(0). 


3. Quantum dot—cavity system 

We now develop a master equation which will al¬ 
low us to derive an analytic expression for <7 hom( t £>) 
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FIG. B.2: Schematic diagram of the system under consid¬ 
eration; a quantum dot with pump-level |p), excited state 
| e), and ground state | g) couples to a cavity mode with 
strength g and detuning A. We consider only the zero pho¬ 
ton and one photon manifolds of the cavity mode, which 
decays with rate k. The \p) —> \e) transition has rate 
a, while the spontaneous emission process of |e) —> \g) 
has background rate F_g. The excited state undergoes 
pure-dephasing with rate 7 . Adiabatic cavity elimination 
results in an effective three-level system with a modified 
spontaneous emission rate T. 


in the presence of time-jitter and pure-dephasing. We 
follow Ref. [H] and model the quantum dot (QD) as a 
three-level-system, with crystal ground state | <?), sin¬ 
gle exciton state |e), and pump level |p), having ener¬ 
gies huig , hui e and hujp respectively. The cavity mode 
is described by creation and annihilation operators c 1 
and c , and has frequency w c . The system is depicted in 
Fig. (B.2). In a rotating frame the QD-cavity system 


is described by the Jaynes-Cummings Hamiltonian 


Hjc = hA\e)(e\ + hg(\e){g\ c + \g){e\ c 1 ), (BIO) 


where A = ( u e — uj g ) — w c is the detuning of the 
QD transition from the cavity mode, and g is the 
QD-cavity coupling strength. Relaxation processes 
are added using the Lindblad formalism [35], and the 
master equation describing the QD-cavity degrees of 
freedom p becomes 

^ = -i[H JC ,p] + (L K (c) + L rB (\g)(e\) 

+T 2 7 (|e)(e|) + L a (\e)(p\)) p (Bll) 

where the Lindblad operators satisfy L^(A)p = 
7 (ApA^ — i{AlH,/ 9 }), with a and k the decay rates 
of the punrp-level and cavity respectively. The back¬ 
ground spontaneous emission rate of the QD is Tb , 
and the rate 7 describes pure-dephasing of the QD 
excited state level. 

In the limit of weak QD-cavity coupling and/or 
strong cavity decay, the cavity can be adiabatically 
eliminated from equations of motion describing our 
system. Formally, we require 7 to t A,T,g with 
7tot = 7 + |(k + r), and provided we consider the 
initial state p( 0 ) = \p)(p\ | 0 )( 0 | with | 0 ) the vacuum 
state of the cavity mode, the dynamics can be well 


approximated by the master equation [28j 
^ t =- i [A\E){E\,p\+ 

[L V (\G){E\) + L 2l {\E){E\) + L a (\E)(P\))p, (B12) 
which is Eq. |l]) in the main text. 


4. Photon Indistinguishability for the 
QD—Cavity System 


We now use Eq. (B12) to calculate the two-photon 


interference probability, Eq. (B 8 ). To proceed, we 


note that in the far field we can make the replace¬ 
ment a(t) —1 <r(f) with er = \G)(E\ [JB] in Eq. (B 8 ). 


Then, to calculate the second order correlation func¬ 


tion G 


( 2 ) 

HBT 


(h,t 2 ) = 


a t (ti)a t (t 2 )a(n)a(Q)) 


^(a^ 1 (ti)^(t 2 )a(ti)a(t 2 )), we make use of quantum 
regression theorem to write [iDl | 


t + r) = —T G® T (t, t + r). (B13) 


For r = 0 we find G^ T (t,t) = 0 since a 2 = (a^) 2 = 
0 , and as such t + r) =0 and we can set 

£/( 2 )(0) = 0 in Eq. (B 8 ). This reflects that for the 
theory presented here we have strictly one (or less) 
excitation in the system at any time. 

We now calculate the two-photon coales¬ 
cence probability expressed in Eq. (B9). To 
begin we consider the uncorrelated probability 
<? ( 2 ) (7,Q) -t (o- t (ti)cr(ti))(CT t (t 2 )o-(Q))- The quan¬ 
tity (a^(t)a(t)) = Tr(p(t)crlcr) is just the e xcited state 
population at time t, and from Eq. (B12|) we have 


(a\t)a(t)) = 0 (f) 


a 


T-a 


(e~ at - e- r * 


) (B14) 


where the Heaviside theta function (0(f) = 0 for t < 0 
and 0 (f) = 1 for f > 0 ) has been introduced to en¬ 
sure no excitations are present before emission events. 
From the quantum regression theorem the first order 
correlation function G^){t,t + r) obeys the equation 
of motion 

^G ( 1 ) (t,f-|-T) = -^ 7 +ir-|-fA^G ( 1 ) (t,t-|-T) (B15) 


with initial condition G ( - 1 )(f, f) = (o’f(f)cr(t)), which 
gives 


G { 1 ) (f,f + r) = (a t (f) t r(t))e- ( 7 + 5 r+lA)|T| . (B16) 


Finally, performing the integrals in Eq. (B8) we arrive 
at Eq. ([ 3 ]) in the main text. 


5. Exciton—phonon coupling 

To explore the influence of phonons seen in our data, 
a weak exciton-phonon coupling time convolutionless 




























master equation technique is used [28] . To second 
order in the exciton-phonon coupling strength, and 
within the Born-Markov approximation, the master 
equation for the complete QD-cavity system (i.e. be¬ 
fore adiabatic elimination) becomes 

^ = -l[H JC ,p\ + (L K (c) + L r (\g)(e\) 

+ L 2 7 (|e)(e|) + L a (\e)(p\)) p + /C p h(p), (B17) 

where the new phonon-induced dissipator is given by 

/»00 

IC P b(p) = - dsTi ph [H I ,[H I (-s),p(t)p ph \], (B18) 
Jo 

where Tr p h denotes a trace over the phonon modes. 
The interaction Hamiltonian is written 

Hi = |e)(e| + b k ) (B19) 

k 

where b' k is a creation operator for a phonon mode 
with wave-vector k , and gk describes its coupling 
strength to the QD exciton. The interaction pic¬ 
ture interaction Hamiltonian is defined by H(—s) = 
e~ lHoS Hie lH ° s , where H 0 = Hjq + H p h, with 
phonon Hamiltonian H p h = u>kb\bk and W& the 

frequency of mode k. Finally, we assume a ther¬ 
mal state for the phonon density operator: p p h = 
e -^ fl p h /Tr p h( e -^ fl ph), with /3 = H/ksT and T the 
sample temperature. 

The strength of the QD-phonon coupling is char¬ 
acterised by the spectral density, defined as J(w) = 


'^2k9k^ u> ~ and which for excitons in QDs has 
been shown to be adequately described by the func¬ 
tion 


J(w) = rjoj 3 4 5 6 7 8 9 10 1 exp [ — (w/w c ) 2 ], (B20) 


where 77 captures the overall strength of the interac¬ 
tion determined by material parameters, and w c is the 
photon cut-off frequency a The behaviour of the 
phonon dissipator in Eq. (B18) in different parameter 


regimes has been discussed in detail elsewhere US¬ 
ED]. The parameters used to obtain improved fits to 
the data in the main text (the solid orange curves in 
Fig. [ 4 ]) are 77 = 0.032 meV -2 and w c = 1.3 meV, 
while the constant pure-dephasing rate was reduced 
to 7 = 1 peV. These parameters correspond to 
phonons contributing approximately 40% of the de¬ 
phasing on-resonance. We note that the other param¬ 
eters in the model were adjusted to 1 /Tb = 730 ps and 
<7 = 34 peV in order that the Tj times as a function 
of detuning were well reproduced. 


The density operator p entering Eq. (B171 contains 


both QD and cavity degrees of freedom. When relat¬ 
ing the field operator a(t) to the QD-cavity system, 
we have a choice to consider QD emission or cavity 
emission, making respectively the replacements a —> c 
or a —> \g){e\ in the field correlation functions. Our 
data was better described by cavity emission, which 
we attribute to the high Purcell factor of our QD- 
cavity system. 
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